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Abstract. - Energy dissipation in a nonequilibrium steady state is studied in driven quantum 
Langevin systems. We study energy dissipation flow to thermal environment, and obtain a general 
formula for the average rate of energy dissipation using an autocorrelation function for the system 
variable. This leads to a general expression of the equality that connects the violation of the 
fluctuation-response relation to the rate of energy dissipation, the classical version of which was 
flrst studied by Harada and Sasa. We also point out that the expression depends on coupling form 
between system and reservoir. 
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Introduction. — Recent developments in nonequilib- 
rium statistical mechanics have clarified fundamental as- 
pects of nonequilibrium fluctuations of work, power flux, 
heat absorbed etc [1-6]. The fluctuation theorem (FT) 
is one of the most remarkable discoveries in nonequilib- 
rium statistical mechanics. This theorem quantifies the 
probability of negative entropy, which can be important 
fbr short measurement times in small systems, and pro- 
vides a precise statement of the second law of thermo- 
dynamics. The relation between the transient version of 
FT and the Jarzynski equality has been demonstrated and 
clarified [5]. Both the FT and the Jarzynski relation [4] 
have been tested experimentally in systems, such as mi- 
Cromechanically manipulated biomoleculcs [7,8], colloids 
in time-dependent laser traps [9,10], and optically driven 
single two- level systems [11]. Studying robust properties 
valid in far-from-equilibrium regime is obviously impor- 
tant for understanding the general structures of nonequi- 
librium statistical mechanics. 

Another important aspect of fluctuations is fluctuation- 
response. In the linear response regime, a fluctuation- 
dissipation relation (FDR) relates a response function with 
an auto correlation function for physical quantities at equi- 
librium [12-14]. However, the FDR is generically violated 
if the system is driven into a nonequilibrium state beyond 
this regime. Several recent studies considered extensions 
of FDR to far-from-equilibrium regime [15-17]. In Refs. 
[15, 16], the violation function was introduced to general- 



ize the Einstein relation, and its validity was experimen- 
tally studied. Harada and Sasa considered the relation- 
ship between the degree of violation of FDR and the rate 
of energy dissipation in over-damped Langevin dynam- 
ics, and found an equality valid in a far-from-equilibrium 
regime [17]. The equality was first derived for the following 



-1x{t)=F{x{t),t)+Ti{t)+ef{t), 



(1) 



where x arc the coordinates, rj is the white Gaussian noise 
satisfying {ri{t)ri{u)) — 2'-fkQTS{t — u), and F{x{t),t) rep- 
resents a force dependent on space and time. Let C{t) be 
the autocorrelation function for the velocity fluctuation 
without perturbation, e = 



Cit) = {[x{t)~v,][x{0)^vs]), 



(2) 



where (...) denotes an ensemble average over thermal 
noise, and Vs is the average velocity of the particle. When 
the external perturbation is finite £ ^ 0, the response of 
the velocity obeys the linear response form [14] 



m 



dt'xit,u)f{u). 



(3) 



Then the following equality was derived in Rcf. [17] 
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[C{u;)-2kBTx'{u;)]}, (4) 
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where / is the rate of energy dissipation. C{uj) is the 
Fourier transformation of C{t) and x' {'^) is the real part 
of the Fourier transformation of xi^)^)- In. an cquiUb- 
rium state with no dissipation flow / = 0, the FDR 
C{ijj) = 2k'QTx{'jj) is satisfied. On the other hand, in 
a nonequihbrium state, the degree of violation of FDR is 
related to the rate of energy dissipation. This equality 
was generalized to correlated thermal noise by Deutsch 
and Narayan [18]. Recently an experimental test of Eq.(4) 
was performed in an optically driven colloidal system [19]. 
In general, it is difficult to conduct direct experimental 
measurements of energy dissipation fiow. Eq. (4) suggests 
that the fiow can be obtained with measurable functions 
in Langcvin dynamics. Thus Eq.(4) is also of practical im- 
portance, providing a new protocol for measuring energy 
dissipation flow in driven Langcvin dynamics. 

In this paper, we consider a driven quantum Langevin 
dynamics and derive a wider class of relations that in- 
cludes the quantum version of Eq. (4) . Quantum Langevin 
dynamics has a wide variety of applications in areas such 
as electronic circuits, superconducting tunnel junctions, 
and electronic systems in semiconductors. [20]. Many dif- 
ferent systems can be mapped onto a simple driven quan- 
tum Langevin equation. We derive a general relation, 
which reproduces the Callen-Welton FDR in an equilib- 
rium state [21]. Our approach provides a unified method 
for investigating energy dissipation flow. This provides a 
new consideration of quantum energy dissipation with re- 
spect to Langevin dynamics, even if it includes nonlinear 
couplings between the system and the reservoirs. 

Quantum Langevin Equation. We consider a 

driven quantum Langevin system described as 



nit) = ^ + v{x,t) 

2m 



^fji 



L 

l2mf 



meujj 



Xi 



\ix V 
me,ujj ) 



(5) 



where {m,x,p\ refer to the system degrees of freedom, 
and {xi,pi,mi^uii\ refers to the reservoir. Those vari- 
ables satisfy the commutation relations [x, xi] = [x(, xt] = 
0, [a;,p] = ifi, and [xi^pgi] = ihdi^i' The potential term 
V{x,t) drives the system into a nonequihbrium steady 
state. Detailed form is not provided here. The coupling 
constant between the system and bath oscillators {A^} is 
switched on at time tini = —00. The initial density matrix 
is assumed to be of the product form pini = ps®PB., where 
S and R refer to the system and the reservoir, respec- 
tively. These matrices are equilibrium distributions. The 
reservoir's density matrix is pn = e~'^^^/Tr[e~'^^^] for 
/3 — l/{kBT). By eliminating the bath's degrees of free- 
dom, we obtain a quantum Langevin equation easily [20], 
which is expressed as 



dt^ 



dVjx.t) 
dx 



duj{u)x(t — u) + rj{t), 



(6) 



where 77 and 7(i) represent a noise term and the mem- 
ory kernels, respectively. These terms control dissipation 
effects from the bath. The properties of the noise and dis- 
sipation are completely determined by the initial condition 
of the bath. We define the spectral function 



A 



ragUJi 



(7) 



Then, the dissipation kernels and noise correlations are 
given by 



2 /-^ , J{uj) 
jit) — — doj coswi. 



{V{t)v{u)) = -/ dc.e— (*-") — J(|c.|)(l + /H) 



where f{uj) = l/(e'"'" — 1). The Fourier transformation 
of the memory kernel is defined as 7(cj) = /„ di7(i)e"^*. 
Then, the real part of the Fourier transformation 7'(w) is 
expressed as 



7'H = JiH)/H 



(8) 



Response Function. — In the Langevin dynamics, 
we consider the relationship between the response func- 
tion and energy dissipation. The response function we 
consider is defined as the response of x against a pertur- 
bation ^ef{t)x. The formal expression of response func- 
tion ^(^7 u) is calculated from the standard linear response 
derivation [14]. In the first order of e, we obtain the devi- 
ation of density matrix from the unperturbed one as 



Spit) 



ih 



dufiu)Uit,tini) [xiu),pitini)]U^it,ti 



(9) 



where pit) is the density matrix at time t, and Uit,tini) 
is the time-evolution operator defined as 



Uitf,t, 



exp^ 



dtnit) 



(10) 



The operator x(w) is U\u,tini)xUiu,tini). We take 
tirii ~> —00, and consider the deviation of the velocity x 
which has the form (3) . Then, we immediately obtain the 
response function, which is expressed using the retarded 
Green function G^^(i,M) as 

xit,u) = -Gl,it,u) = ^eit-u){[xit),<u)]). (11) 

Here, (...) denotes an average over the initial state. 

Energy Dissipation and the Callen-Welton FDR. 

— We consider an energy dissipation fiow from the system 
into the reservoir. By taking the derivative of the system's 
Hamiltonian with respect to time, we get the heat current 
operator 



Z = - V —pxe + V 



A 



1 



2miujf m 



ixp + px), (12) 
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driving field. In order to simplify it further, we consider 
the following quantity for the Green function 



Fig. 1: Schwinger-Keldysh contour. 



1 fT/2 i-r/2 

G^B(tj)= lim - / / dtdt'G'XB{t,t')e"^'e-"^\ 



where positive current flows from the system into the 
reservoir. To arrive at the average current at the steady 
state, we employ the technique of Keldysh green function. 
We use conventional notations for the Green functions, 
which are defined for arbitrary operators, A and B, as 

G\s(.t,t') = -'-{{A{t)B{t')) + {B{t')A{t))), 

GTBitX) = -Je(±(i-t'))([S(i),A(t')]), 

where G^5(t, f), G^i3(i,<')> and G\g{t,t') are the 
Keldysh, retarded, and advanced Green function, respec- 
tively. We set the initial time to be tint = —t/2. Using 
the green function, an average current is calculated as 



lim — 

r — >oo T 



r/2 



(I'W) = -?E^^7|-Gl..(ii,i2^ 



r/2 

_d_ 

dti 



dt{I'it)} + Io, 



(13) 
(14) 



where /o = liuir^^ ^ J^'^^^^dt ■^{x'^{t)) J^ duj"/'{uj)/Tr, 
which is the contribution from the second term in 
Eq.(12). We assume that the initial density matrix of 
the system is an equilibrium distribution for the Hamil- 
tonian without the potential part V{x,t), i.e, p'^/2m + 
(x^/2) J2e ^£ / {"iT^e^j) ■ This set up for initial state enables 
us to use the Wick theorem to compute the Green func- 
tions. Perturbation expansions for contour-ordered Green 
functions are performed along the Schwinger-Keldysh con- 
tour depicted in Figure 1. Using the Langreth rule for 
making the Keldysh Green function in Eq.(14) [22], we 
readily derive the average current 

/ = /o + lim — / / dtdt' 

X {GL(t,t')S^(i',i)+G'L(i,0S'^(i',t)}-(15) 

In Eq.(15), the function S'"' and S"^ are self-energy terms 
from the reservoirs, which calculated from the free green 
functions for reservoir's Hamiltonian. These are written 



I]fc(i',t) = ^^ / duj-f'{uj) '^^^^^— cos{iu{t' - t)) , 

^ ^ 7T Jo ^anh(^) ^ ^ '^ 

-2 (9 f°° 

j:^{t\t) = — Q{t-t')— duJj'{uj)cOs{uj{t' -t)) . 
n dt' Jq 

The formula (15) is valid for an arbitrary time-dependent 



where a represents ?■, a, and k. When translational in- 
variance in time is satisfied, this reduces to the usual ex- 
pression of the Fourier transformation. Let x'('^) denote 
— Re[Gj3.(ti')]. Then, a straightforward modification leads 
to 



/ 



1 



dw^' {lo) 



2t:v15{uj) + C{lo) 



X'{uj)nuj 



tanh(^) 



(16) 



where we used the integral by parts to remove the term 
/o. Here, G(w) = {iri/2)G'l,j^{uj) - 2ttv15{uj), and Vs is the 
average velocity. In an equilibrium state where no net 
energy dissipation flow exists, the Callen-Welton FDR, 
C{ijj) — x'(^)^w/ tanh(/3?itj/2) is satisfied [21]. In the 
limit ?i — > and by inserting ^'{io) = 7, the above equa- 
tion reproduces Eq.(4). 

Concluding Remarks. — Several general principles 
exist for nonequilibrium phenomena. In the linear re- 
sponse regime, the validity of the Onsager reciprocity and 
the Green-Kubo relations have been established. The fluc- 
tuation theorems and Jarzynski equality are being inves- 
tigated in numerous models and experiments, and they 
seem to be the exact relations valid arbitrarily far from 
equihbrium. Eqs.(4) and (16) are also interesting, show- 
ing that the degree of FDR violation is related to energy 
dissipation in Langevin dynamics. This would be impor- 
tant, because Langevin dynamics is ubiquitous in realis- 
tic systems. Quantum Langevin dynamics is believed to 
have wide applicability in many realistic systems, includ- 
ing metallic tunnel junctions with capacitances and su- 
perconducting junctions [20,23-25]. In those cases, the 
system and reservoir variables represent the variables in 
electrical circuits, and the external force is realized by an 
electric current. Circuit-realization of a driven harmonic 
trapped particle was proposed in Ref. [6]. The quantum 
case for this would be an interesting relevant system. 

In the present study, Eq.(15) is the key equahty for 
deriving the result (16). To derive Eq.(15), linear cou- 
pling between the system and thermal environment was 
critical. Although linear couplings should be the dom- 
inant contributions in most realistic systems, nonlinear 
couplings, where a nonlinear function of x couples with 
the bath variables, are also possible [26]. It is possible to 
generalize Eq.(15) to nonlinear coupling cases. In general, 
when we change the coupling form in the Hamiltonian (5) 
as \ixix — > Xexef{x), a different expression of average 
current from Eq.(16) is obtained as 



1 
I = - 



dujj '{lo) 



~2^ff 



H + 



TiwRcfG'; 



■,M] 



tanh(iif^ 



(17) 
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This means that the expression depends on types of cou- 
pUng form. It would be important to figure out how Eqs. 
(4) and (16) are generahzed, if we consider other types of 
reservoirs and dynamics. Energy dissipation of spin dy- 
namics would be an important problem to be studied. 

Another intriguing problem might be on higher order 
fluctuations of energy dissipation flow. Nonequilibrium 
fluctuations increases in time unlike equilibrium ones. To 
study characteristics of fluctuations, it is convenient to use 
the technique of counting statistics [27,28], which is equiv- 
alent to the protocol for obtaining fluctuation theorem in 
quantum systems [29-31] and work distribution [32]. It is 
possible to derive an explicit form of characteristic func- 
tion, which generates not only average current but also 
any orders of cumulants of dissipation flow. It can re- 
produce the present result (16). Systematic derivatives of 
the characteristic function with respect to a counting field 
generates any orders of fluctuations. We hope that this 
study encourages further studies on energy dissipation at 
far-from-equilibrium conditions in quantum systems. 
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